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Abstract 

Excited-state relaxations in molecules are responsible for a red shift of the 
absorption peak with respect to the emission peak (Franck-Condon shift). 
The magnitude of this shift in semiconductor quantum dots is still unknown. 
Here we report first-principle calculations of excited-state relaxations in small 
(diameter < 2.2 nm) Si nanocrystals, showing that the Franck-Condon shift 
is surprisingly large (~ 60 meV for a 2.2 nm-diameter nanocrystal) . The 
physical mechanism of the excited-state relaxations changes abruptly around 
~ 1 nanometer in size, providing a clear demarcation between "molecules" 
and "nanocrystals". 
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The Stokes shift commonly observed in molecules QXJ] and ionic solids 0, has its origin in 
excited-state atomic relaxations. When an electron-hole pair is created by optical excitation, 
the final state has approximately the same atomic configuration as the initial state (Franck- 
Condon principle). Prior to emission, however, the system can relax to a new configuration 
with lower total energy. Recombination occurs from the relaxed atomic configuration, lead- 
ing to a red-shift of the emission lines with respect to the absorption lines (Franck- Condon 
shift). The Franck-Condon shift in molecules and ionic solids can be as large as several eV. 
The magnitude of excited-state relaxations in semiconductor quantum dots, on the other 
hand, is still controversial. Quantum dots grown by colloidal chemistry methods range in 
size between between 10 2 and 10 4 atoms, so excited-state relaxations could be significant. 
Continuum models based on the effective-mass approximation, however, have predicted a 
Franck-Condon shift of only a few meV |3|,[| . 

A schematic diagram of the relevant electronic energy levels as a function of the generic 
configuration variable R is shown in Fig. 1. The minimum-energy atomic configuration of 
the quantum dot in the electronic ground state (GS) is different from the minimum-energy 
atomic configuration of the dot in the singlet (S) or triplet (T) excited states. At low 
temperature and in the absence of light, the quantum dot is in the ground-state geometry 
Res'- The lowest-energy allowed optical transition takes the system into the singlet excited 
state, as the lower-energy triplet state is optically inactive. Since the excitation of an 
electron-hole pair occurs on a much faster time scale compared to atomic vibrations, the 
atomic configuration immediately after absorption is still the ground state configuration 
Ros, as indicated by the vertical line in Fig. 1. Absorption can actually occur into several 
vibrational states associated with the singlet electronic state, leading to a broadening of 
the absorption line. The exciton relaxation then proceeds according to a few characteristic 
times: (i) The spin-flip time Tfi ip , which is the time required for the exciton to switch from 
the singlet state to the triplet state, (ii) the recombination times r^'J (with rf ec <C t^ c ), 
which include radiative and non-radiative recombination paths, and (iii) the relaxation times 

S T 

r rl ' x in the singlet and triplet states, which are the characteristic times for the dissipation of 
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the vibrational energy. The relaxation times r r i x are ultimately determined by the coupling 
of the quantum dot with the environment. If T rl '^ -C rj ec (and Tfu p < tJ £C ), the quantum dot 
relaxes to the lowest-energy triplet excited-state configuration (Rxs) before the electron and 
the hole can recombine. In this case the total Stokes shift is AE = AEpc + AE$t, where 
AEpc is the Franck-Condon shift and AE$t is the exciton exchange splitting. Different 
vibrational states may be involved in the emission process, leading to a broadening of the 
emission peak fT]j2|]. 

In this Letter we investigate, using ab-initio density-functional methods, the excited- 
state dynamics of Si quantum dots, and calculate the ensuing Franck-Condon shift. We 
find that for small Si quantum dots (1.0 — 2.2 nm diameter) the Franck-Condon shift is 
surprisingly large. For example, in the case of a ~ 2.2 nm-diameter quantum dot we predict 
a Franck-Condon shift of ~ 60 meV, versus an electron-hole exchange splitting of only 
~ 8 meV. By analyzing the physical mechanism responsible for the Franck-Condon shift 
we are able to identify two physically distinct regimes: for sub-nanometer Si clusters the 
dominant mechanism is the stretching of a single Si-Si bond upon electronic excitation, while 
for larger Si nanocrystals the Franck-Condon shift originates from a change in the overall 
shape of the nanocrystal in the presence of an electron-hole pair. This distinction provides 
a clear demarcation line between the molecular regime and the nanocrystal regime. 

The calculations were performed using ab-initio density-functional theory in the local 
spin density (LSD) approximation. We used ultra-soft pseudopotentials to describe the 
electron-ion interaction, and the plane-wave representation (with an energy cutoff of 150 
eV) to describe the Kohn-Sham orbitals. Triplet excited states can be calculated within 
density-functional theory by minimizing the total energy of the system in the triplet spin 
configuration. We find that in practice this approach works very well. For example in the 
case of the silane molecule SiH 4 , we calculate a triplet excitation energy of 8.1 eV, compared 
to the value of 8.7 eV obtained using diffusion Quantum Monte Carlo ||. In the bulk limit, 
we expect our approximation for the excited-state energy to converge to the LSD band gap 
of Si, which is ~ 0.68 eV lower than the experimental band gap. So we expect the calculated 
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excitation energies to differ from the experimental excitation energies approximately by a 
constant energy shift 0. 

Our approach to the calculation of excited states is applicable only to the lowest-energy 
triplet state (see Fig. 1), where the electron excited to the LUMO and the electron remaining 
in the HOMO have parallel spin, i.e. \T) = | if). Calculations of the excited-state singlet 
energy surface would require the handling of a two-determinantal wave function, which is 
beyond the reach of simple density-functional theories (see however Ref. [0). For comparison, 
we have also performed excited-state calculations on the "mixed" energy surface given by 
\M) = | T I)- ^ can be shown that, for a given atomic configuration, the energy difference 
E M — E T is one half of the singlet-triplet splitting E s — E T (see Ref. 0). 

We consider here nearly spherical Si nanocrystals centered on a Si atom. The initial 
atomic configuration (before atomic relaxations) is obtained by cutting out a sphere from 
a bulk Si crystal. The Si-Si bond length is taken as the bulk LSD bond length (2.33 A). 
The surface atoms having three dangling bonds are removed, while the remaining surface 
dangling bonds are passivated by H atoms. The nanocrystals considered here range in size 
from 29 Si atoms (10.3 A diameter) to 275 Si atoms (21.7 A diameter). 

The calculation of the Franck-Condon shift is carried out in 4 steps: (i) First, the ground- 
state atomic configuration is obtained by minimizing the total energy of the nanocrystal with 
respect to the atomic positions, as dictated by quantum-mechanical forces. This step gives 
the ground-state total energy E gs (Rgs)- (ii) Then we excite an electron-hole pair in the 
triplet state, and calculate the excited-state energy in the ground-state geometry: E t (Rgs)- 
The difference E T (R GS ) — E GS (TL GS ) is the triplet excitation energy E exc (not to be confused 
with the absorption energy E a b s shown in Fig. 1). (iii) Next, we relax the atomic positions 
on the triplet excited-state energy surface, thus obtaining the excited-state total energy 
in the excited-state atomic configuration (XS): E T (TL X s)- (iv) Finally, we calculate the 
ground-state total energy in the excited-state geometry: E GS (Rxs)- The energy difference 
E T (Rxs) ~ E GS (Rxs) is the emission energy E emi . The Franck-Condon shift AEpc is then 
given by 
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AE 



FC 



E t (Rgs) — E GS (~R GS ) — E T (R X s) — E GS (R X s) = E exc — E emi . (1) 



AEpc can be further decomposed into an excited-state contribution AExs = E t (Rgs) — 
E T (R X s) and a ground-state contribution AE GS = E GS (R X s) - E GS (R X s) (see Fig. 1). 
Note that AExs and AEqs are total-energy differences between different atomic configura- 
tions of the same system. Thus, we expect the accuracy of the calculated Franck-Condon 
shift to be comparable with the accuracy of calculated vibrational energies (a few % error 
in the case of bulk Si). 

The results for Si nanocrystals are summarized in Table I. Note that, in the size range 
considered here, the ground-state contribution to the Franck-Condon shift AEqs is larger 
than the excited-state contribution AExs, particularly for the smaller nanocrystals. This 
difference reflects the reduced curvature and increased non-parabolicity of the excited state 
energy surface compared to the ground-state energy surface. Martin et al. [|3| proposed 
a simple model, based on the envelope-function approximation and empirical deformation 
potentials, to estimate the Franck-Condon shift in Si nanocrystal. They predicted that for 
a nanocrystal with an excitonic gap of 2.3 eV (corresponding to the largest nanocrystal 
considered here, see Table I) the Franck-Condon shift would range from 9 to 21 meV, 
depending on the parameters of the model. Using a similar continuum model, Takagahara 
et al. [§]] obtained a Stokes shift of ~ 7 meV for the same nanocrystal. Our ab-initio 
calculations show that these models significantly underestimate the Franck-Condon shift of 
Si quantum dots. 

The Franck-Condon shift is very large in the small Si29 H36 cluster, where the distortions 
due to the electronic excitation are large. Figure 2(a) shows the Si-Si bond-length distribu- 
tion n(L) of this cluster, both in the ground-state geometry and in the relaxed excited-state 
(triplet) geometry. We see that n(L) is similar in the ground state and in the excited state, 
except for a single Si-Si bond that is stretched by about 15% in the excited-state geometry. 
The HOMO and LUMO wave functions in the excited-state geometry are strongly localized 
around the stretched Si-Si bond, as shown in the insets in Fig. 2(a), and the correspond- 
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ing energy levels are well inside the band gap of the Si nanocrystal, thus determining the 
large Franck-Condon shift of 2.9 eV. Using tight-binding total energy calculations Allan 
et al. |8| predicted that in small hydrogen-passivated Si nanocrystals excitons can become 
self-trapped in a (meta) stable state localized at the surface. Our ab-initio calculations show 
that excited-state relaxations lead to the spontaneous formation of a stretched bond in the 
interior of the nanocrystal, and therefore this effect should depend weakly on the type of 
surface passivation. Hirao || calculated the Stokes shift of Si29 H 36 nanocrystals, finding a 
value of 0.22 eV, over an order of magnitude smaller than our result. This difference may 
be due to the lower energy cutoff used in the plane- wave expansion of Ref. 0. 

Figure 2(b) shows the bond- length distribution of the Si§7 H_7 6 nanocrystal. The distri- 
bution is centered around the bulk Si-Si bond length, with very little differences between 
the ground state and the excited state. In fact, we find that the Franck-Condon shift in 
this nanocrystal (and in larger nanocrystals) is due to a change in the overall shape of the 
nanocrystal, from spherical to ellipsoidal, upon electronic excitation. This change of shape 
leads to a splitting of the states at the top of the valence band (which are degenerate in the 
Trf representation) and thus to the Franck Condon shift. The insets in Fig. 2(b) show that 
the electron and hole wave functions in the excited-state configuration are delocalized over 
the entire nanocrystal. 

As shown in Fig. 1, the exchange contribution to the Stokes shift is given (at low 
temperatures) by the singlet-triplet splitting AE S t evaluated at the ground-state geometry 
{Res'}- We therefore calculate the exchange contribution as: 



AE 



ST 



E M (R GS )-E T (R GS ) . (2) 



The results are shown in the last column of Table I. We find that the electron-hole exchange 
splitting is significantly smaller than the Franck-Condon shift, even for the largest nanocrys- 
tal considered here. The singlet-triplet splitting of Si quantum dots had been calculated in 
the past using a variety of empirical approaches, including the effective-mass approximation 
ill , the tight-binding approximation H, and the semi-empirical pseudopotential method 
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11 1 . Our ab-initio calculations are in good agreement with the pseudopotential calculations 



ofRef. 11 



The Stokes shift of Si nanocrystals has been measured using both optical and thermal 
methods p~0j , [12|JT3|1 . It is difficult to compare our results directly with experimental data, 
because Si nanocrystals synthesized by colloidal chemistry techniques are typically larger 
than those considered here, and because the absorption coefficient at the absorption edge is 
small. Using selective laser excitation Calcott et al. |10| and Kovalev et al. [p~|] measured the 
red shift of the photoluminescence onset with respect to the excitation energy. For example, 
for an excitation energy of 2.41 eV, Calcott et al. [1(J measured a red shift of 23 meV. 
The photoluminescence onset was attributed to zero-phonon emission, and the observed 
Stokes shift was interpreted as due to the electron-hole exchange splitting. Our calculations 
suggest that in small Si nanocrystals the Franck-Condon shift becomes progressively more 
important, until it dominates over the exchange splitting. 

In conclusion, we have shown by excited-state density-functional calculations that the 
Franck-Condon shift in small Si nanocrystals is larger than previously thought. We have 
found that in sub-nanometer clusters the Franck-Condon shift originates from the stretching 
of a Si-Si bond, while in larger nanocrystals it is due to a change in the overall shape of the 
nanocrystal upon electron-hole excitation. 

This work was supported in part by the DOE Computational Materials Science Network 
grant DE-FG02-02ER45972, NSF grant DMR9803768, the US DOE under contract DE- 
AC05-00OR22725 with the Oak Ridge National Laboratory, managed by UT-Battelle, LLC, 
and the William A. and Nancy F. McMinn Endowment at Vanderbilt University. 
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TABLES 

TABLE I. Calculated Franck-Condon shift AEpc and electron-hole exchange splitting AEst 
(in eV) of a few hydrogen-passivated Si nanocrystals. Also shown are the ground-state and ex- 
cited-state contributions to the Franck-Condon shift: AEpc = AExs + AEqs- 



Nanocrystal 


Diameter (A) 


AE XS 


AE GS 


AE FC 


A Est 


Si29 H36 


10.3 


0.79 


2.13 


2.92 


0.051 


Si 8 7 H 76 


14.8 


0.12 


0.20 


0.32 


0.021 


Sii47 H100 


17.6 


0.09 


0.12 


0.21 


0.014 


Si275 Hl72 


21.7 


0.03 


0.03 


0.06 


0.008 
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FIGURES 

FIG. 1. Schematic diagram of the ground-state (GS) and excited-state singlet (S) and triplet 
(T) energy surfaces of a semiconductor quantum dot. Light is absorbed (E a b s ) by exciting the 
quantum dot from the ground state to the optically active singlet state. Emission occurs (E em i) 
from the optically inactive triplet state, leading to a resonant Stokes shift of the emission line. The 
Franck-Condon contribution to the Stokes shift is the sum of AExs and AEgs- 

FIG. 2. Bond-length distribution in the ground-state geometry (dashed lines) and in the triplet 
excited-state geometry (solid lines) for the Si2gH36 nanocrystal and the Sis7H76 nanocrystal. The 
insets show the electron (LUMO) and hole (HOMO) single-particle wave functions in the ex- 
cited-state configuration. The arrows denote the calculated bulk Si-Si bondlength. 
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